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Studies of the hydrogen bond. III 
Molecular orbital treatment of some saturated and unsaturated ethers 
By Inca FiscHEerR 


Introduction 


The concept of the hydrogen bond is generally accepted among chemists, 
although the precise nature of hydrogen bonding is still not fully understood. 
Some experimental evidence on the hydrogen bond, and viewpoints as to its 
character were reported in two preceeding papers of this series (HHRENBERG 
and Fiscuer [1]), in which the influence of adjacent double bonds on the 
ability of a “heteroatom” (N, O) to be “acceptor” in hydrogen bond formation 
was studied. In particular, the unsaturated compound divinyl ether, and the 
corresponding saturated compound, diethyl ether, were compared in respect of 
their admixture with water. Several experimental facts suggest that diethyl 
ether forms hydrogen bonds relatively easily, whereas divinyl ether has no 
tendency to enter into hydrogen bond formation. 

In this paper some results of theoretical calculations of hyperconjugation 
energies, charge distributions, and bond-orders in these ether molecules are 
given, and some conclusions regarding hydrogen bonding by ethers are reached. 
The calculations were made using the molecular orbital (m.o) method (for a 
rewiew of this method: see CouLson [2]). 


Method 


The usual approximation made in applying the m.o. method to a conjugated 
system is to consider only the interactions of the z-electrons, and to neglect 
the o-electrons entirely. We have followed this convention in the calculations 
of this paper. For a molecule of » atoms, J, w..., and n z-electrons the j:th 
m.o., expressed in the LCAO (linear combination of atomic orbitals) approxima- 
tion, 1s 

v= Sted 
u 
where g, is a normalised, hydrogenlike 2p orbital of the atom mw, and aj, 1s 


a parameter to be determined by a variational method. The normalisation 
condition of the m.o. is 


foi vide = Daj a Sin =], 


A, 
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where S;,=J yi g.dt is the overlap integral. The energy of the electron, 
E;, 1s = i| 
E;= Jyi Hyde 


where H is the one-electron Hamiltonian appropriate to a self-consistent field. 
The variational treatment gives in the usual way the determinantal equation 


|Hi,—ES8i,|=0 (1) 


where Hi, = {gyi Hgdt. Hj,= a is a Coulomb integral and Hi, = Bin an 
exchange integral. Equ. (1) is of the v:th degree in #. To each root, Hj, » 
corresponds a set of parameters, 7,, which determine the 7:th m.o. In the 
ground level two electrons occupy each of the m. o:s corresponding to the lowest 
energy values. 

When overlap is included the electron charge on the atom A, q,, and the 
x-bond order between the atoms A and wm, pin, are defined by Coulson (private | 
communication) as: 

ag = Dei Yi () = occupied m. o:s) 


I 


\ 
Yiu = » Sip Lji- 


A 


where 


Now it is convenient (see e.g. WHELAND [3] and MUuLLIKEN et al. [4]) to in-— 
troduce a quantity y: 


Yiu = Biw — $ Sin (as + Ou) (2) 
y is connected in a simple way with the bond energy, as shown below. Let 
ae — K=z; Oy Oe = 20 (3) 


where 6, is a measure of the difference in electronegativity between the atom 
jw and the carbon atom. Introducing (2) and (3) in equ. (1), gives: 


re aE 07 Viw a Siu (2 a oF aia Ou) a 
Yiu aia 7m (Z ae 0) ae On) z2+2 On Soe neice 0. (4) 


i ee ee eae RE em Tha boride CG Mice Oo UR Min 4 


oY fej, 10)” eh es, “ea Fe eee ae 1S fe Set Se wl ae-e.) 4 tej een a Pits ere 


KEqu. (4) may be solved if we know the numerical values of y and 6. Neither 
of these quantities are known exactly, but we can obtain probable values of 
y and 6 from their connexion with experimentally known bond energies, letting 
6 attain values reasonably compatible with the different proposed electronega- 
tivity scales (see e.g. Pautine [5], Mutriken [6]). This connexion is found by 
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considering a m.o. extended over two atoms only, one of which is carbon. 
The highest root of (4) is 


il 
db = =a eS 2 
&@ (ins 8?) Ly, S, Von Sait Yul 


where the index C is omitted. If we assume that a, and a, have the same 
values for bonding and non-bonding electrons, the binding energy, D’, is given 
by the expression: 


, 2 2 2 
D! = e+ O, —2(ag—2) = LY Se ey. (eee) |} 


D’ is approximately equal to the energy difference between a double and a 
single bond between the atoms in question, provided that all the quantities 
are referred to the same interatomic distance (for a detailed discussion on this 
point, see RoBEerts and SKINNER [7]). This is obtained through the construc- 
tion of the Morse curves D = Do (1—e74""~"0')? of the bonds in question. Do 


is the dissociation energy, and the constant a is calculated from a = Vk/2 Do, 
where & is the force constant. The values of these parameters used in the 


_ present paper are given in Table I. 


Table I 
D r+ 108 k-10° | 
Bond type to) co} 

oP keal cm dyn.cm 

ONE acre sy Stee | 145 1.33 9.9 

eC row ene aisusnsonekeisys 83 1.54 3.8 

(Oh (0) eee niko 156 OAL 12.0 

(Gea EN eee oon ca ethene 80 1.43 | 4.4 

OR C7 i] os Doe ae a eee 98 1.09 4.9 
Oa Sen eee 110 0.96 Shi ial 


The overlap integrals were evaluated using the 2p, atomic wave functions 
given by Duncanson and Coutson [8]. The following screening constants were 
used: carbon 3.12, oxygen 4.70, hydrogen 1.16. As z-orbitals of the H3- and 
H,-groups we used the linear combination of atomic orbitals given by MULLIKEN 
et al. [4]. As to the values of 6 and y, we assigned reasonable values to 6 and 
calculated the corresponding values of y from equ. (5). In order to survey 
the influence of a variation in 6 on charge distribution and bond order we 
performed the calculations with different sets of 6, y-values. These are sum- 
marised in Table II. The values used for the CC-bonds are: 


fos 1.354 S = 0.290 y = 2.132 eV 
%o = 1.54 » S = 0.213 y=1.147 » 
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Table II 
| CO CH OH 
| = 
| |» = 143 | pee Lae r = 109 p= 10. hr =e 
| S= 0.138|} S= 0.141 S'= 0.586 S= 0.339 | S= 0.149 
—20,eV | yeV yeV +20, eV yeV —20op eV yeV yeV 
0.5 | 1.220 1.268 0.2 3.370 1.2 321 — 
0.8 1.167 L223 0.3 3.367 is) 3.076 0.744 
1.0 1.128 1.180 0.5 3.000 Le 3.041 0.591 
ee IS a RO0 1.124 0.7 3.343 
125 0.956 1.013 : 
Results 


Dimethyl ether. The structural formula of this molecule may be written 
1.09 1.43 


H,=C—O—C=H;. We consider hyperconjugation involving two 2p, electrons 
from the oxygen and one pair of z-electrons from each of the H,;=C triple 
bonds. In this as well as in the two following cases the secular equation may 
be factorised by considering the symmetry of the molecule, which permits a 
sub-division into orbitals symmetrical, y*, and antisymmetrical, y%, in the nuclei. 
The hyperconjugation energy, R, 1s given by: 


ao (Coe eee aia Oop 


For this simple molecule the calculations were performed with several com- | 
binations of dy and odo values. Table III contains the computed values of | 
hyperconjugation energy, charge distributions, and bond orders. We _ have | 
tabulated the net charges, Q., 7.e. the sum of z-electronic and nuclear charges 
(the contribution from the o-electrons is not taken into account). It is seen | 
that for a given 06, the charge on the oxygen atom and the bond order feo | 
are almost unaffected by variations in dy. The hyperconjugation energy, how- 

ever, varies considerably with both 6, and dy. Furthermore, Table III shows | 
that when dy = 0.3 eV we obtain a larger, formal negative charge on H than | 
on C. As this is rather improbable, we limit ourselves in the following to | 
calculations with dg > 0.3 eV. 


Diethyl ether. This molecule, which has the formula 


1.09 1.54 1.48 
Hj=C=0—0—0=C =H. 
l I 


2 2 


was treated analogously to dimethyl ether. The part of the hyperconjugation : 
energy due to the inclusion of the oxygen atom, is given by the expression: | 


5 
R=2 D4 + 400—4 (24 + 24). 


j=1 
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Table III 
Dimethyl ether 


rr 
S il 
Ae | 
2 Q | 
H H Q | 
eV NX ‘ Qo PHC Peo keal | 
S 
26, =— 0.5 eV 
0.2 —0.053 —0,022 +0.150 0.961 0.273 
022 LO, 2 8.09 
0.3 —0.039 — 0.036 40.148 0.962 0.270 7.72 
0.5 — 0.010 —0.062 +0.146 0.963 0.264 6.98 | 
25,=— 0.8 eV | 
0.2 —0.047 —~ 0.020 0.134 0.966 0.258 6.26 | 
0.3 — 0.032 — 0.034 +.0.132 0.966 0.255 5.91 
0.5 —0.004 —0.061 0.130 0.967 0.250 5.22 
0.7 +0.024 — 0.088 +0.128 0.967 0.244 4.53 
| 20,=— 1.0 eV 
0.2 — 0.042 —0.019 +0.122 0.969 0.246 4.92 
03 — 0.028 — 0.033 +0.121 0.969 0.243 4.60 
0.5 +0.001 — 0.060 +0.118 0.970 0.238 3.95 
0.7 +. 0.029 —~ 0.088 +0.116 0.969 0.233 3.30 
26,=—12eV 
02 | =0,036 —0.019 +0.109 0.972 0.232 3.53 
0.3 — 0.022 —0.032 +0.108 0.972 0.230 3.23 
0.5 | +0.007 — 0.060 0.106 0.973 0.225 2.64 
0.7 40.035 —0.087 0.104 0.973 0.220 2.03 
| 2 do =— 1.5 eV 
0.2 — 0,026 —0.018 0.087 0.978 0.208 41 
0.3 —0.011 — 0.032 0.086 0.978 0.205 1.16 
0.5 40.017 — 0.060 +4.0.084 0.978 0.201 0.65 
0.7 +.0.046 — 0.088 +.0.083 0.978 0.196 0.14 


The numerical results obtained for diethyl ether are given in Table IV. For 
this molecule the values of the parameters dy and 6) were varied within a 
smaller range than in the preceding case, partly because the calculations are 
much more cumbersome than those for dimethyl ether, partly because it was 
found, that the quantities, which are of special interest for the present in- 
vestigation, 7.€. Qo, Peo, and R, are almost identical with the corresponding 
values of dimethyl ether. 


Divinyl ether. For this molecule, 


1.35 1.42 ; 
aL aa ie ; He ; 
| 
H H 


the hyperconjugation energy is obtained from the expression: 
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3 
R=2 dia — 428 t+ 460. 


j=1 


The computed values are summarised in Table V. Some computations were} 
made with a 5 % lover value of Seo in order to see how much this variaviGg 
changes the results. The influence was found to be negligible. 


Table IV 
Diethyl ether 


Qn’ Qe On Q Qo PH'O' Pole PHC Peo 


| 267 =0.50eV 26,=—10eV _f 

+0.071 —0.070 —0.002 —0.057 +0.116 0.982 Osa 0.955 0.234 3.91 
20y=0.7eV 26,=—10eV 

+0.100 —0.099 +0.026 —0.083 -+0.113 0.980 0.183 0.953 0.228 3.28 9 
264=0.50eV 26,=—15eV : 

+0.071 —0.071 +0.015 —0.056 +0.082 0.982 0.178 0.963 0.197 0.66 | 


20, = O07 eV 26, == ToeVv | 
| +0.100 —0.099 +0.042 —0.083 +0.080 0.979 0.184 0.961 0.192 0.17 


Table V 
Divinyl ether 


Qo 

| eV MG 
| Soo = 0.141 
Lea Que —0.181 +0.015 + 0.332 0.906 0.415 14.76 
ES — 0.166 +0.022 + 0.287 0.918 0.387 11.40 

1.0 0.153 + 0.026 + 0.255 0.927 0.366 9.08 
lee Slee —0.140 +-0.029 + 0.222 0.936 0.343 6.75 
peariae —0.117 + 0.031 +0.171 0.951 0.3038 3.32 8 
| Seo = 0.135 
0:5 — 0.176 +0.015 +0.327 0.907 0.410 14.94 | | 
Pe, Peetip lien +0.030 +0.168 0.952 0.299 3.59 1 


Ethyl vinyl ether. We thought it of some interest to compute bond orders: 
and charge distribution of the asymmetrical ethyl vinyl ether. The calculations: 
were made with the following parameter values: 26) =— 1.5 eV, 26y = 0.5 eV, I 
de — ac’ = 0.3 eV; C belongs to the vinyl group and C’ to the ethyl group.. 
a fg values of net charges and bond orders are indicated in the: 
ormula: 
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—0.123 +0.032 +0.140 —0.095 —0.118 +0.113 
, 


H,C——CH—_0—C’_(’ == 


0.951 0.303 0.222 | 0.184 0.972 


0.954 He, 
40.046 


It is seen that they agree well with the corresponding results for the symme- 
trical ethers. 


Either-water complexes. As is shown by Counson and FiscHer [9] the simple 
m.o. method is unreliable when the interatomic distance is increased to more 
than 1.675. In most hydrogen-bonded molecular complexes the distance be- 
tween the hydrogen of the one molecule and the oxygen of the other is near 
to this critical value. Though well aware of this fact we have made a ten- 
tative investigation of the hydrogen bond using the m.o. method in the crude 
form described above. Although one should not overemphasize the value of 
the results, it may be assumed that they give a correct qualitative description 
of the hydrogen bond. 


\ 
We consider the system O'..... H—O, where O’ belongs to an ether. We 
1.6 1.0 


assume that the three atoms lie on a straight line with the OH interatomic 
distances given in the formula. “‘Molecular’’ orbitals obtained from a combina- 
tion of the atomic orbitals which have rotational symmetry around the nuclear 
axis lead to the charge distribution and bond orders collected in Table VI. It 


Table VI 


Ether-water complex 


LG : 9 | 
goo Son Qo’ Qn Qo Po'R PHO 
eV eV 
5 145) + 0.050 + 0.232 — 0.282 0.256 0.943 
» ee + 0.041 +0.188 — 0.229 0.247 0.954 
laa ihe7/ +0.033 + 0.267 — 0.300 0.218 0.945 
» 1.5 + 0.032 + 0.242 — 0.273 (Ore LIE 0.952 
» 1) +0.030 +0.191 — 0,221 0.200 0.965 


is seen from the Table that the o-bond order between the hydrogen of the 
water and the oxygen of the ether, 2.e. the bond order of the hydrogen bond, 
is around 0.2 for all the ao- and a-values considered. A purely electrostatic 
interpretation of the hydrogen bond should require this bond order to be naught. 


Discussion 


Firstly, some conclusions may be drawn from the numerical results, with 
regard to the m.o. method. The value chosen for 0, affects to some extent 
the charge on the atom mw and the order of the bonds to this atom. The in- 


501 


I. FISCHER, Studies of the hydrogen bond. III 


fluence on other charges and bonds decreases rapidly with the distance. The») 
influence from removed parts of the molecule is so small that a methyl group ) 
and an ethyl group affect the ether oxygen to the same extent (Tables IIT | 
and IV). The similarity between the electron distribution of the asymmetrical! | 
ether and the symmetrical ones also emphasizes the fact, that removed parts) 
of the molecule do not seriously influence one another. This result shows, that. 
in the treatment of homologous molecules it is justified to use the same valueny 
of the parameters 6 and y. It also justifies the use of diatomic m. o:s in the ) 
calculation of corresponding 6- and y-values (equ. (5)). The y-values from (5}) 
do not agree, however, with MuLirKken’s' suggestion that y/S = constant, , 
which may depend on the approximations involved in both relations. The facts, , 
well known from experiments, that bond energies and bond dipole moments are} 
almost constant within a homologous series are in good agreement with the: 
narrow range found for the influence of a variation of 6,. — The effect on. 
bond order of a variation of 6 and y has earlier been discussed by CouLson’ 
[10], who varied 6 and y independently. 

The value of the hyperconjugation energy is much more sensitive to a change: 
in a than are charge and bond order. MUuLiikEN [11] has shown that the «! 
value of carbon varies from molecule to molecule, not only with a change from 
triple to double bond but also with the radicals bound to the atom. This 
emphasizes the extreme difficulty in choosing the “best” a value. When, as 
in the present case, the effects of hyperconjugation are studied it therefore : 
seems better to draw conclusions from the calculations on bond order and. 
charge distribution than from those of delocalisation energy. 

The results obtained are subject to many uncertainties. We shall now discuss, 
the influence of some of these factors. A questionable point is the expressions | 
for the o- and z-orbitals of the Hs- and H,-groups. We have adopted the: 
usual expressions though the overlap integral of the CH z-bond turns out te 
be twice as large as the overlap of the corresponding o-bond. This result (also 
obtained by CovuLson [12]) suggests that the z-bond should be stronger than | 
the o-bond, and needs further investigation. The point is, however, not of } 
especial importance to the present results, as was shown by performing the | 
computations with a 20% lower value of Don. This change in D’ did not |} 
affect the net charge on the oxygen atom at all and only caused a small | 
change in the charges on the carbon and hydrogen atoms. Nor are errors of | 
a few percent in the assumed bond energies of importance. | 

The values of the overlap integrals depend on the screening constants, c, for | 
the atomic wave functions, on the internuclear distances, 7, and on the steric | 
configuration of the molecules. It was found that the use of a c-value ap- | 
propriate to O* instead of O (in the main computations we used a c-value be- | 
tween these) diminished the net charge by 4 %. The use of different screening 
constants for saturated and unsaturated ethers would not, therefore, diminish | 
the ratio between their Qo-values by more than 1 %. Even less is the effect | 
of an error in r of the order 0.02 A. The steric configuration of divinyl ether | 
has been studied using the electron diffraction method by Bastransen [13], to 
whom I am much indebted for performing that investigation. According to | 
him, there seems to be a deviation from the planar configuration in divinyl] 
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ether, the C=C bonds making an angle of 15—20° with the plane of the C—O 
bonds. This will not, however, diminish the S-value by more than 5 %. This 
variation in S is considered in Table V. 

In the present investigation we have not accounted for orbital hybridisation. 

According to BasTIANSEN [13] the O—C—O angle is approximately the same 
in diethyl and divinyl ether. This means the same degree of hybridisation of 
the oxygen 2p orbitals occurs in the two compounds. 
In general, the m.o. method will give comparable results in the case of 
similar molecules. This is certainly true in case of dimethyl ether and diethyl 
ether. As for diviny! ether there is reason to believe that this molecule is 
rather different from the saturated ethers. In our opinion, the main difference 
lies in a difference in the values of a> and ae and, consequently, of 65. From 
a study of the a:-values computed by MuturKen [11] the conclusion can be 
drawn that the difference of —265 between the saturated and unsaturated 
ethers most probably is less than 0.5 eV, divinyl ether having the higher value. 
a i computation for ethyl vinyl ether the difference was assumed to be 
ove. 

The values given in Tables III, IV, and V show clearly that the net charge 
and the bond order of the oxygen atom of saturated ethers are less than the 
corresponding quantities of unsaturated ethers. The hyperconjugation energy 
shows the same trend, though the result is less conclusive on this point. 
Analysing the influence of the various uncertainties considered above, we find 
that they will lower the ratio between the Qo- and eco-values of divinyl ether 
relatively to those of diethyl ether by less than 3 %. Consequently, we may 
draw the following conclusion: an adjacent double bond changes both the charge 
on an etheral oxygen atom and the C—O bond order, 7%.e¢. the polarisability 
of the oxygen atom. 

In order to compare some of the theoretical values with experimental ones 
we choose as a probable case the one with —26, of divinyl ether equal to 
1.2 eV and the corresponding value of diethyl ether 1.0 eV. The tabulated 
QYo-value of divinyl ether is + 0.22 and that of diethyl ether + 0.11. Because 
of the C—O o-bonds the oxygen atom is negatively charged in both ethers 
though the computed z-electron distribution will diminish this charge, more in 
case of divinyl than of diethyl ether, the difference corresponding to a ~ 0.2 
Debye smaller dipole moment. This agrees well with the experimental differ- 
ence [1]: 1.23—1.06=0.17 Debye. The computed peo-values, 0.34 of divinyl 
and 0.23 of diethyl ether, means a bond length difference of 0.02 A in good 
agreement with electron diffraction values. 

We may now apply our results in respect of the nature of the hydrogen 
bond. The electrostatic energy of the complex considered above will be roughly 
proportional to Qo’ as there is no reason to assume different O—O distances 
in the water complexes of the two ethers, which have almost the same sterical 
configuration around the oxygen atoms. Using the values of Qo from Tables IV 
and V we find a difference in electrostatic energy between the complexes with 
water of diethyl and divinyl ether, which for all the — 20 -values considered 
is less than 2 kcal/mole. This difference is less than half the energy of the 
hydrogen bond, which is known to vary between 4 and 8 kcal/mole. Thus, a 
purely electrostatic interpretation of the hydrogen bond is incomplete. We also 
have to consider polarisation effects. The formation of the hydrogen bond will 
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disturb the hyperconjugation with the rest of the molecule of the 2 pz electrons. 
of the ether oxygen. As co is 1.5 times as large in divinyl ether as in diethyl | 
ether the activation energy for the water complex formation should be con- 
siderably larger in the case of divinyl ether. Thus, divinyl ether cannot enter 
into hydrogen bond formation because the oxygen atom is too little polarizable. 
This explains why the solubility of water in divinyl ether is as low as in| 
benzene [1]. Diethyl ether, however, can form hydrogen bonds relatively easily as 
the CO z-bond order of this ether has about the same value as the O’...H 
o-bond order in the hydrogen bond, see Tables IV and VI. | 

We are led to conclude that the energy of the hydrogen bond depends not 
only on electrostatic attraction but also on a certain delocalisation, 7. ¢. po- 
larisation, of the 2 electrons of the “hydrogen acceptor” atom. Quantitative 
conclusions about the degree to which these two effects contribute to the 
energy of the hydrogen bond cannot be drawn from the present investiga- 
tion, because the values obtained for the delocalisation energies are too un- 
certain. 


This work was started under the direction of Professor C. A. CouLsoNn 
London. It is a pleasure for me to express my sincere thanks for his stim- 
ulating guidance during my stay in London and for his continued, helpfu! 
interest in my work. My thanks are also due to Fil. kand. Mrs. B. Mosack, 
who assisted with the numerical calculations. 


SUMMARY 


Charge distribution and bond orders are calculated with the simple m. o. 
method for the following ethers: dimethyl, diethyl, divinyl, ethyl vinyl. The. 
influence of the values chosen for the Coulomb integrals is discussed. 

The m.o. method is also applied to the hydrogen bond of a molecular 
complex of ether-water. It is found that the 2 p electrons of the etheral oxygen 
atom are delocalised to some extent through the hydrogen bond formation. 
Thus, the energy of the hydrogen bond is not purely electrostatic in origin. 
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